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The Algol triple system spatially resolved at optical wavelengths 
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Interacting binaries typically have separations in the milli-arcsecond regime 
and hence it has been challenging to resolve them at any wavelength. However, 
^ I recent advances in optical interferometry have improved our ability to discern 

r^ I the components in these systems and have now enabled the direct determination 

Y^', of physical parameters. We used the Navy Prototype Optical Interferometer to 

2 ' produce for the first time images resolving all three components in the well-known 
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Algol triple system. Specifically, we have separated the tertiary component from 
the binary and simultaneously resolved the eclipsing binary pair, which represents 
the nearest and brightest eclipsing binary in the sky. We present revised orbital 
elements for the triple system, and we have rectified the 180-degree ambiguity in 
the position angle of Algol C. Our directly determined magnitude differences and 
masses for this triple star system are consistent with earlier light curve modeling 
results. 



Subject headings: astrometry 
— stars: individual(Algol) 



binaries: eclipsing — techniques: interferometric 



1. Introduction 

Algol (/? Per, HR 936), the prototype for a well known class of eclipsing binaries, has 
well over 200 years of publishe d observations ava ilable for study. The eclipsing nature of the 
system was first suggested by iGoodrickd (1 17831 ) when he stated the light variations could 
result from "... the interposition of a large body revolving aroun d Algol ..." . A long susp ected 
third component in Algol was spectros copically confirmed by IStruve fc Sahadd (119571 ) and 
Ebbighausenl (119581 ). ISoderhjelml (Il980l ) provides a review of Algol which is a useful starting 
point for summarizing this triple system. We will refer to the B-type primary star in Algol 
as Algol A, t he K type secon dary star as Algol B, and the more distant Am companion as 



Algol C (see iRicharda Il993l. for th e exact spectral types). It is thus a hierarchical triple 



system as de fined bvlEvana (119681) . I Wade fc Hjellmind (119721 ) made the first radio detection 
of Algol and iLestrade et al.l (119931 ) determined that the radio emission in Algol comes from 
Algol B using multi-epoch Very Long Baseline Interferometry observations. Algol was one 
of 12 radio stars used to link the HI PPARCOS optical refer ence frame to the International 
Celestial Reference System (ICRS) (IKovalevskv et al.lll997). T he 1.86 year orbit of Algol 
C was resolved b y speckle interferometry (jLabeyrie et al.lll974[ l and optical interferometry 
( jPan et al.l 119931 ) but these investigations could not resolve the close binary. Also, these 
results suffer from a 180° ambigu ity in the absolute position angle (PA) due to a lack of 
phase information (JLabeyridll970l ). The Fourth Catalog of Interferometric Measurements of 
Binary Stars has a more complete listing of these observation^!]. 

This considerable body of knowledge still leaves room for additional exploration, and 
permitted inconsistent descriptions of the orbital elements of the triple system. It would 



^http://www. usno.navy.mil/USNO/astrometry/ optical-IR-prod/wds/int4 (JHartkopf et al.ll2001al) 
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appear that the orbital elements listed in ISoderhjelnj (jl980[ ) and iPan et al.l ( 119931 ) would 
adequately describe the Algol triple system. But these elements include a 180° difference 
in the position angl e of the ascending nod e when compared with the orbital solution in the 
Hipparcos catalog ( jLindegren et al.l 119971 ). Algol's role as one of a small number of radio 
stars used to link the HIPPARCOS optical reference frame to the ICRF requires resolution of 
this inconsistency. Algol can also serve as a position angle calibrator for optical and near-IR 
interferometers and this strengthens the case for a resolution of this disagreement. 



Recently, Algol was observed in the near IR (K^ = 2.133 /im) wit h the CHARA array an d 
at 5 GHz with the European VLBI Network during December 2006 (jCsizmadia et al.l 120091 ). 
The CHARA array with approximately 200 m baselines has a similar spatial resolution to 
the NPOI at optical wavelengths with 64 m baselines. The CHARA data resolve the close 
binary but ICsizmadia et al.l (|2009[ ) make no mention of detecting Algol C. Csizmadia et al. 
use their VLBI and CHARA array observations to p roduce an orbit of Al gol A— B with 
an opposite sense of rotation from that determined by iLestrade et al.l (119931 ). We are then 
presented with another inconsistency in published results for the Algol system. 

During 2006 October and November we collected observations of Algol usingthe Navy 
Prototype Optical Interferometer (NPOI), the Very Long Baseline Array (VLBA)g, and the 
Lowell Observatory 42" Hall telescope equipped with the solar— stellar spectrograph. We 
defer discussion of the radio and spectroscopic results for a future paper. The primary aim 
of this project was to resolve the close pair in Algol with the NPOI, and perform an absolute 
astrometric registration of the optical NPOI images to the ICRF phase-referenced VLBA 
images. As the NPOI records visibility squared and closure phase data we can determine the 
position angle calibrati on without a 180° amb iguity and resolv e the inconsistency between 
the orbital elements in ISoderhjelml ( 1l980l ) and iPan et al.l ( 119931 ) and the orbital elements in 
the HIPPARCOS cat alog. We can al s o add ress t he inconsistency i n the direction of the close 
binary orbit between ILestrade et al.l ( 119931 ) and I Csizmadia et al.l (120091 ) . 



In this paper, we report on our NPOI observations which extend the direct knowledge 
of the Algol triple system to optical wavelengths. In §2 we discuss our observations, with 
particular attention to the calibration of the NPOI absolute position angles. In §3 we present 
the astrometric orbits of the Algol A— B and AB— C systems. We conclude with §4 and a 
comparison with light curve solutions and a discussion of the astrometric results. 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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2. Observations and Calibration 

Algol was observed with the NPOI from 1997 October to 2008 Octo ber (see Table El). 



The N POI is a six element optical interferometer, described in detail in lArmstrong et al. 
( 119981 ) ■ De tails regarding t he ge neral NPOI observationa l setup and data recording can 



be found in iHummel et al.l ( 120031 ) and iBenson et al.l (|2003[ ). We combined the astrometric 
siderostats on the center (AC), east (AE), and west (AW) stations with the E6 and W7 
imaging siderostats. The addition of the latter allowed projected baselines of up to 64 m 
in length, separated by about 60 degrees in position angle. Up to three baselines were 
recorded on each of the two spectrometers. We switched between two different four-station 
configurations half-way through the night. Post-processing of the data was performed using 
C. Hummel's OYSTER software package. The calibrator star for the 2006 observations was 
e Persei (HR 1220, V=2.89, B0.5V, parallax (tt) = 6.06 mas), loc ated 9.6° from Algol and 



an estimated diameter (d p,st) of 0.43 mas based on its i? — / color (JMozurkewich et al.lll991 



White fc Feiermanlll987| ). The 2008 observations included e Per, check binary 7 Per (HR 
915), and an additional calibrator for 7 Per, e Cassiopeiae (HR 542, V=3.34, B3111, tt = 7.38 
mas, dest = 0.43 mas). The uncertainties on the estimated diameters are 10%. Since 7 Per 
is 16.2° from e Per, Algol and 7 Per could share the same calibrator, e Cas then serves as a 
secondary check on the calibration. 

To verify the absolute NPOI PA calibration using our 2008 Oct 27 observation of 7 
Per we fit a model to the obs erved squared visibilities and triple p hases in OYSTER and 
imaged 7 Per using DIFMAP (jShepherdlll997l : IShepherd et al.lll994| ). The expected position 
of 7 Per was 163.0 mas at 244.6° with an expected R band magnitude difference of 1.5—1.6 
dPrieur et al]l2003h . We observed 7 Per at a position of 160.92 ± 0.28 mas at 244.95 ± 0.22° 
and a magnitude difference (800nm) of 1.50 ± 0.08. The same result was obtained for both 
choices (e Per or e Cas) of the calibrator stars which verifies the absolute PA calibration of 
the NPOI. 



3. Analysis and results 

The combined visibility data for each night allow the determination of the relative posi- 
tions of the components. Note that the orbital motion of the A— B pair is significant during 
each night's observation. Reliable estimates of the magnitude differences were available in 
the literature and were used as an initial guess in our model fits. The dominant feature of 
a finite magnitude difference between the C component and the combined light of the A— B 
pair in the data is the pronounced sinusoidal variation of the squared visibilities. Super- 
imposed on this variation is a subtle modulation due to the (larger) magnitude difference 
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of the close binary itself. The magnitude differences in the V [AV = 2.92) and Cousins-I 
(A / = 2.63) bands b etween the eclipsing pair and component C was determined previously 
by iPan et al.l (Il993[ ) using the Mark III Stellar Interferometer. An initial estimate for the 
magnitude diff e rence between A and B was determined from the light curve analysis by 
Richards et al.l (119881 ) in the V band {AV = 3.92), and we estimated the val ue for the NPOI 
800 nm filter (A/ = 2.6) using the effe ctive tempera tures and log g values of lRichardd (119931 ) 
and the Kurucz model atmospheres ( jKuruczl 119791 ). As a c heck, the same p roced ure cor- 
rectly reproduces the J-band magnitude difference given by [Richards et al.l (119881). Since 
the diameters are only barely resolved we adopted the values given in JRichards et al.l (Il988[ ) 
for components A, B, and C, converted to angular diamet e rs usi ng the parallax of Algol of 
35.1 mas (distance 28.5 pc, distance modulus 2.27, IESA I (Il997l )). The stellar parameters 
initially used to model the Algol triple system are listed in Table [T] and were kept fixed for 
the fits of the relative component positions to refine the orbital elements. After these initial 
refinements of the orbital elements we also fit for the masses. 

Initial astrometric results (separation and position angle) were obtained by using an 
image for each night to provide an initial guess for the separation and position angle. Algol 
and 7 Per were analyze d in the same manner. I mages of Algo l were made using AIPS 
dvan Moorsel et al.lEggeh . DIFMAP and BSMEM dBuscherl Eggj ) . Figure E] illustrates the 
motion of Algol B over two epochs and the location of Algol C with uniformly weighted 
images made with DIFMAP. This guess was refined in a fit to the visibility data directly. 
As typically done in the reduction of NPOI data, we used a fraction of the CLEAN beam 
(20% in this case) to provide a more realistic estimate of the uncertainty ellipse since the 
formal errors from the fit to the visibility data usually underestimate the true uncertainty in 
the results. These initial astrometric results were used to derive an initial fit for the orbital 
elements. 



For the A— B orbit we fixe d the eccen tricity e to zero ( ISoderhjelml Il980l ). We used 



the photometric light elements 



Kim 



19891) f or the A— B orbit, and the inclination from 



the light curve analysis ( iRichards et al.l 119881 ). A fit was then performed for the position 
an gle of the ascending n ode f l, and semi-major axis a. The Na D lines of Algol B detected 
by iTomkin &: LambertI ( 119781 ) and their orbital elements do verify that we have correctly 
identified the quadrant of the ascending node of Algol A— B. 



For fitting the AB— C orbital elements we began with the orbital elements of iPan et al. 



(119931 ) and then corrected the position angle of the ascendin g node. We solved for all the 



orbit al elements. Again, comparison with spectroscopic data (lEbbighausenlll958l : iHill et al. 
197ll ) verify that the quadrant of the ascending node is correctly identified. 



In the final step, these orbital elements were refined by comparing them with the vis- 
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ibility data for each scan. Our final astrometric results appear in Table [2l Differential 
corrections due to the orbital motion of the A— B pair are included in Table [2] based on 
the orbital elements listed in Table |3i Thes e differential correctio ns were first applied to 
NPOI data for close binaries as described in iHummel et al.l (Il995[ ). The relative positions 
in Table |2] between composite components in the hierarchical triple, i.e. AB-C, refer to the 
photocenter of AB. 

Using the Algol C orbital elements we again solved for the magnitude differences (B and 
C relative to A) and the three elements (a, f2 & P) of the A— B p air. A final fit for th e stella r 
masses was perf orme d using the radi al velocities contained in iTomkin fc LambertI ( 119781 ): 
Hill et all f ll97lh and IhIU et aP f ll993h for Algol A, B, and C. The masses in Table [D served 
as initial guesses for the fits. The results and uncertainties are given in Table Ei Figure [2] 
shows the orbits of the Algol system with the astrometric data. 



Discussion 



Our discussion of the analysis of the interferometric observations will be restricted to the 
new insights we have gained on this well studied system. Our NPOI observations mark the 
first resolution of the Algol system into three components. ICsizmadia et al.l ( 120091 ) reported 
the resolution of the close binary with the CHARA array at near infrared wavelengths but 
without absolute phase calibration. Our o rbital solution i s fully consistent with the pio- 
neering radio interferometric observations of iLestrade et al.l (jl993[ ). We have unambiguously 
determined that the close pair orbit is retrograde, and nearly orthogonal to the plane of the 
wide orbit. The relative angle is 96° ± 5° according to 



cos (j) = cos Zi cos 12 + sin ii sin ^2 cos( Afi) , 



where ii and ^2 are the inclinations of the two orbits, and AQ is the difference between the 
two ascending node angles. This improvement t o the orbital plane o rientation is relevant 



to dynamical studies of hierarchical triples (e.g. JKiseleva et al.lll998l ). This orientation of 
the wide orbit also removes the discrepancy e xpected be t ween the photocenter motion one 
would compute using th e AB— C elem ents of JPan et al.l ( 119931 ) and that found from HIP- 
PARCOS or the orbit of iHeintzl (Il994j). Our retrograde orie ntation of the close binary orbit 
co ntrasts with the progra de orbit of ICsizmadia et al.l (120091 ). The position angle calibration 
of ICsizmadia et al.l (120091 ) depended on VLBI observations made during a radio flare of Algol 
which m ay have complicated th eir analysis due to the changing radio morphology of Algol 



4.1 of ICsizmadia et al.l 120091 ). As the NPOI observations are calibrated to produce an 
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absolute position angle as shown in §2 (with 7 Per as the position angle calibrator) we are 
confident that the retrograde o rbit o f the close binary pair is correct. The 15 GHz VLBI 
observations of iPeterson et al.l ( 120101 ) also agree with the previously determined retrograde 
orbit. 

The absolute position ang le calibrat i on of the NPOI enables a revision of the orien- 
tation of the Algol C orbit of iPan et al.l (119931 ) . Ou r results place the maximum AB— C 
separation in the same quadrar it as found by Heintz (JHeintzl Il994j : iGatewood et al.l Il995l ) . 
The Hipparcos o r bital solution ( JESA Ill997l . Double and Multiple Systems Annex) references 
Gatewood et al.l ( 119951 ) for the quadrants of the longitude of periastron u and the position 
angle of the ascending node Vt, and our observations confirm the accuracy of the Hippar- 
cos orbital elements. This confirmation is important because it links results from both the 
optical and radio reference frames. A 180° reversal of the position angle of the ascending 
node of the AB — C orbit presented he re would create a time variable systematic offset of 
the photocenter (IGatewood et al.lll995l ) that could not be reconciled with the HIPPARCOS 
observations. 

Our interferometric observations have resulted in the first resolved images of the triple 
system and the first directly measured magnitude differences for the three stars in Algol. 
Previous estimates of the magnitude difference for the close binary were ra ade by modeling 



photometric and spectroscopic data (e.g. iKimI Il989l : iRichards et al.l Il988l . and references 



therein) . The V band magnitude d ifferences predicted from these m odels span slightly more 
than on e magnitude: 3 72 ± 0.10 Jwilson et all 1972 lZ^7 ± 0.31 ( ISoderhielmlll98ot ). 3.92 
± 0.88 jRichards et al.lll988h . and 2.71 ± 0.15 (lKimlll989h . Our directly measured V band 
magnitude differences (Table |3]) favor magnitude differences of less than 3. Our AB— C 
magnitud e difference is in excellent agreement with the previously determined value using the 
Mark HI (iPan et al.lll993r) and it is also consistent with an early speckle interferometry result 



( jLabeyrie et al. 



I974J ). We extrapolated the magnitude difference to the center wavelength 



of the Hipparcos Hp filter using the stellar atmosphere parameters from Table [H and used 
the masses given in that table to determine a 18.4 mas amplitude for the motion of the 
photo center. This result is consistent with the Hipparcos orbital solution of 19.0 ± 0.6 mas 
(ESA 1997). The dynamical parallax determined from our full fit to the astrometric data 
and pubhshed radial velocities is 34.7 ±0.6, consistent with the Hipparcos value of 35.1 ±0.9 



mas. 



Our magnitude differences add a directly measured constraint to the results obtained 
from modehng the photometric light curve and spectroscopic data. It may be useful to re- 
examine the modeling of the close binary using magnitude differences derived directly from 
the interferometer measurements. Other bright double and multiple stellar systems will yield 



similar constraints for use with spectroscopic and photometric data. 

The Navy Prototype Optical Interferometer is a joint project of the Naval Research 
Laboratory and the US Naval Observatory, in cooperation with Lowell Observatory, and is 
funded by the Office of Naval Research and the Oceanographer of the Navy. The authors 
would like to thank Jim Benson and the NPOI observational support staff whose efforts made 
this project possible. S. Bosken & G. Shelton of the USNO Library were very helpful in our 
literature search. RTZ is thankful for the support of the Fizeau Exchange Visitor Program 
through the European Interferometry Initiative (EII) and OPTICON (an EU funded frame- 
work program, contract number RII3-CT-2004-001566). This work was partially supported 
by National Science Foundation grant AST-0908440 to Richards. The literature search made 
use of the NASA ADS abstract service, the SIMBAD database maintained at CDS, Stras- 
bourg, France and the JSTOR archives. We thank an anonymous referee for numerous 
suggestions which improved the paper. 

Facilities: NPOI () 
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Table 1. Adopted and Derived System Parameters 



Parameters 


Stellar Component 




A 


B 


C 


Diameter (mas) 


0.77 


0.93 


0.37 


Mass (Mq)'^ 


3.7± 0.3 


0.81± 0.05 


1.6± 0.1 


Teff(K)- 


13000 


4500 


7500 


log g 


4.0 


3.5 


4.5 



Note. 



Values for diameters calculated as de- 



scribed in §3. The log(g) values are from iRichards 
(1l993[ l rounded to match the atmosphere models we 
used. Diameter, T^e and log g were fixed during the 
model-fitting. The masses here were initial estimates 
in the model fitting and the final mass results and 
uncertainties are listed in Table |3l 



Richards! fll993[ ) 



Table 2. NPOI Relative Astrometric Results 







AB-C 




A.-B 


Error Ellipse 


a 




UT Date 


JY 


P 


e 


P 


e 


^maj 


^tniti 


</- 


* 




(yrs) 


(mas) 


(deg) 


(mas) 


(deg) 


(mas) 


(mas) 


(deg) 




1997 Oct 16 


1997.7901 


100.63 


309.22 






0.996 


0.334 


2.0 




1997 Oct 17 


1997.7928 


100.94 


309.96 






0.790 


0.366 


168.0 




1999 Mar 4 


1999.1700 


12.11 


228.86 






1.360 


0.396 


133.9 




2006 Oct 19 


2006.7976 


50.56 


300.60 


2.08 


226.49 


1.006 


0.256 


9.6 


0.78 


2006 Oct 20 


2006.8004 


51.56 


299.26 


1.98 


48.10 


1.056 


0.270 


5.0 


0.14 


2006 Oct 23 


2006.8086 


53.38 


299.58 


2.54 


41.01 


1.002 


0.248 


181.5 


0.18 


2006 Oct 27 


2006.8195 


55.30 


301.52 


1.25 


248.71 


0.572 


0.256 


-1.4 


0.57 


2006 Oct 28 


2006.8223 


56.32 


301.79 


1.29 


216.29 


0.660 


0.222 


169.6 


0.93 


2006 Oct 29 


2006.8250 


57.36 


301.01 


2.05 


54.25 


0.460 


0.286 


153.2 


0.27 


2006 Oct 30 


2006.8277 


57.38 


302.45 


1.95 


230.14 


0.668 


0.232 


173.4 


0.62 


2006 Oct 31 


2006.8305 


58.26 


302.82 






0.440 


0.272 


146.0 




2006 Nov 1 


2006.8332 


58.78 


301.49 


1.80 


43.12 


1.014 


0.228 


123.5 


0.32 


2006 Nov 2 


2006.8360 


59.21 


302.73 


2.06 


227.45 


0.944 


0.232 


119.1 


0.67 


2006 Nov 3 


2006.8387 


60.16 


302.22 






0.924 


0.204 


129.3 




2006 Nov 4 


2006.8414 


60.55 


301.92 


1.52 


40.09 


0.886 


0.220 


125.4 


0.36 


2006 Nov 5 


2006.8442 


60.95 


303.31 


2.33 


227.33 


1.506 


0.238 


100.0 


0.72 


2006 Nov 6 


2006.8469 


62.08 


302.41 


1.43 


75.06 


0.866 


0.242 


106.1 


0.06 


2008 Oct 27 


2008.8208 


83.00 


306.52 






0.564 


0.292 


11.1 





Note. — Col. 1: UT Date of observation; Col. 2: Julian year at 7 hours UT; Cols. 3-6: separation, position angle (from 
north through east) for the AB-C and A-B components, respectively; Col. 7: Semi-major axis of error ellipse; Col. 8: Semi- 
minor axis of error ellipse; Col. 9: Pos ition angle of error ellipse; Col. 10: Close binary orbital phase $ using light elements 
2441773.49 + 2.8673285 * E JKimll 19891) . * of 0.0 = primary eclipse, 0.5 = secondary eclipse. 

''The error ellipse is the uncertainty in the location of the position vector. For component C this is with respect to the AB 
photocenter. 
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Table 3. Obital Solution And Component Parameters 



Orbital 


A-B 


A-B 


AB-C 


AB-C 


Element 


Soderhielm (1980) 


This work 


Pan et al. (1993) 


This work 


a (mas) 


2.2 ± 0.1 


2.3 ± 0.1 


94.61 ± 0.22 


93.8 ± 0.2 


i (deg) 


81.4 ± 0.2^ 


98.6^ 


83.98 ± 0.09 


83.7 ± 0.1 


n (deg) 


132 ±4 


47.4 ± 5.2 


312.26 ± 0.13 


132.7 ± 0.1 


e 








0.225 ± 0.005 


0.225 ± 0.005 


oj (deg) 






310.29 ± 0.08 


310.8 ± 0.1 


To (JY) 




1973.2471*^ 


1987.3689 


1987.3689 


To (JD) 




2441773.491^ 


2446931.4 ± 1.5 


2446931.6 ± 0.1 


P (days) 


2.8673 


2.867328 


680.05 ± 0.06 


679.85 ± 0.04 


P (years) 






1.8619 ± 0.0002 


1.8613 ± 0.0001 


TTdyn (mas) 








34.7 ± 0.6 


Magnitude Differences 








Components 


Am(550nm) 


Am(800nm) 






A-B 


2.70 ± 0.3 


2.20 ± 0.3 






A-C 


2.8 ± 0.2 


2.6 ± 0.2 






AB-C 


2.9 ± 0.1 


2.7 ± 0.1 






Masses (Mq) 










M(A) 


3.7 ± 0.2 








M(B) 


0.8 ± 0.1 








M(C) 


1.5 ± 0.1 










ssl) 








^Richards et al. (IS 





^Minimum light of primary eclipse 



^i > 90° used to indicate retrograde motion as defined by iHeinta (119781 ). 
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Fig. 1. — An image of the Algol triple system made from the NPOI data of 2006 Oct 29. Algol 
C is the component in the upper right hand corner. The inset shows a close-up image made 
from the NPOI observation of 2006 Oct 30 and emphasizes the motion of Algol B between 
the two epochs. To guide the eye the approximate positions of Algol B are indicated at each 
epoch by a filled black circle. The uniformly weighted restoring beam is shown in the lower 
right hand corner. 
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Fig. 2. — The large figure illustrates the AB— C orbit. A vector from the origin indicates 
the periastron point. The inset figure shows the A— B orbit with a vector from the origin 
indicating the position of primary eclipse minimura light. The astrometric results of Table [2] 
are plotted with the astrometry of iPan et al.l (119931 ) rotated by 180° and the orbital elements 
in TableO Uncertainty ellipses are 20% of the CLEAN beam for the NPOI data. Astrometric 
positions are fit to the individual 30 second scans. The plotted positions for A— B are 
computed at UT07:00 on the date of observation between 2006 Oct 19 and 2006 Nov 06. 
Arrows show the direction of the orbital motion on the sky. 



